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ABSTRACT 


Convolution  integral  oquations  have  a  simple  representation  in  the 
form  of  a  Laplace  transform,  and  so  solution  follows  by  transform  inversion. 
Usually,  the  solution  is  expressed  as  an  infinite  integral  which  can  only  be 
evaluated  by  numerical  integration.  Often,  the  tramsform  involves  the  ratio 
of  fianctions  of  hypergeometric  type.  In  piievioQBi.,papayc  (sae  ’’Thu.  Piidfe' Table 


Hi  11  III  mily* approximations  to  an  inverse  Laplace  transform^vhen  the  trans¬ 
form,  which  involves  the  exponential  integral,  is  approximated  by  rational 
functions.  The  approximate  solution  is  a  sum  of  exponential  functions,  and 
numerics  are  presented  to  show  the  efficiency  of  our  technique. 
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INTRODUCTIOM 


Numerous  problems  in  applied  mathematics  require  the  solution  of 
integral  equations  of  the  convolution  type.  It  is  well  known  that  equations 
of  this  kind  have  a  rather  simple  representation  in  the  form  of  a  Laplace 
transform  so  that  the  solution  depends  upon  finding  an  inverse  Laplace  trans¬ 
form.  Now  it  may  be  that  the  Laplace  transform  involves  functions  for  which 
rational  approximations  are  available,  and  it  is  of  interest  to  study  the 
nature  of  approximations  to  inverse  Laplace  transforms  when  the  transforms 
are  approximated  by  rational  functions. 

In  previous  studies  [l,2j,  we  have  been  concerned  with  rational 
approximations  to  higher  transcendental  functions.  The  approximations  are, 
of  course,  useful  to  evaluate  the  transcendental  functions  as  in  the  sense  of 
table  making.  However,  a  more  cogent  reason  for  investigating  such  approxi¬ 
mations  is  that  they  should  be  useful,  for  example,  to  invert  Laplace  trans¬ 
forms,  and  in  general  to  further  simplify  the  evaluation  of  mathematical  solu¬ 
tions  to  applied  problems  already  expressed  in  closed  form. 

In  two  recent  studies  (j5,43,  we  have  investigated  problems  in  ap¬ 
plied  mechanics  and  aerodynamics  whose  solution  requires  inversion  of  a 
laplace  transform  which  involves  ratios  of  functions  of  hypergeometric  type. 

In  the  cases  studied,  the  hypergeometric  functions  were  the  modified  Bessel 
functions  of  the  second  kind  Ky(z)  .  With  the  aid  of  results  developed  in 
[2j,  we  derived  economic  approximations  for  the  inverse  functions.  In  the 
present  paper,  we  present  a  similar  study  where  the  transcendent  involved  is 
the  exponential  integral. 


I.  SOLUTION  OF  A  CONVOLUTION  INTEGRAL  EQUATION 

For  our  particular  study,  we  take  an  example  frcan  a  paper  by 
Friedlander  [s]  which  requires  the  solution  of 


z 

y*  k(z-t)g(t)dt  =  g(z)+f(z)  , 

o 

k(z)  =  2f(z)  =  -2(2+2)-2  .  (l.l) 
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Using  the  theory  of  Laplace  transforms,  we  have 


G(P)  =  ' 

» 

F(P)  =  -  i  -pEl(-2p)  ,  -Ei(-2p)  -  f  t-^e-'*'dt 

2p 

0(p)  .  I  }  , 

2f2pe^Zi(-2p) 


(1.2) 


Here  F(p)  ,  G(p)  and  K(p)  are  Laplace  transforms  of  f(z)  ,  g(z)  and 
k(z)  ,  respectively.  Also  -El(-z)  is  the  usiial  notation  for  the  exponential 
integral,  ^y  -the  inversion  formula  for  Laplace  transforms. 


6(z) 


1 

2ni 


c+i® 


/ 

C-i“ 


eP^(p)dp 


(1.3) 


where  c  >  0  and  c  lies  to  the  right  of  all  singularities  of  G(p)  . 

In  the  remainder  of  this  section,  we  develop  some  results  which  lead 
to  an  integral  representation  of  g(z)  .  The  results  obtained  are  also  useful 
in  connection  with  our  process  for  approximating  g(z)  .  We  first  show  that 
G(p)  has  no  poles  in  the  conplex  plane  cut  along  the  negative  real  eixis,  and 
that  G(p)  is  analytic  every\diere  except  on  this  axis  for  at  the  origin  it 
has  a  branch  point.  To  prove  this,  we  first  note  that 


z  “  k  k 

-Ei(-z)+(Y+ln  z)  =  r  t"l(l-e''‘^)dt  =  -  ^  ,  (l.4) 

0  k=l 


X 

Ei(x)  =  £  f^e'^dt  ,  X  >  0  ,  (l.5) 
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Ei(z)  =  Ei  ('-ze-^^^+iTT  =  I  Ei  (-ze^^  ^  +Ei  ,  (l.6) 


-ze^Ei(-z)  '-i  2Eo(l^li-z"^)  =  ^  (-)^k;z"^  ,  jz|-^co 

k=o 


jarg  z|  <  3tt/2  . 

(1.7) 


In  (1.4)  ;  Y  is  Euler's  constant,  and  in  (1.5)  the  integral  is  a  Cauchy 
principal  value.  Also  in  (l.7)  and  in  the  sequel,  standard  generalized  hyper- 
geometric  notation  is  used.  Thus 


Xa2^,a2, . .  .Sp  X 


IT  (‘i),-** 


■  Z  -  ,  (.)v  =  . 


k=o  ^ 

TT  (»i)|^k: 

i=l 


r(a) 


(I.e) 


For  further  discussion  of  hypergeometric  functions,  see 
Nov  let 


h(z)  =  1  +  i  (2z)e22Ei(-2z) 


(1.9) 


and  consider  the  contour  described  by  the  figure  below. 
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I 
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Here  let  be  the  small  circle  with  radius  r  ,  and  T2  the  large  circle 

with  radius  R  .  The  complex  plane  is  cut  along  the  negative  real  axis  so 
that  along  BC  ,  arg  z  is  n  ,  while  along  EF  ,  arg  z  is  -n  .  Then  the 
number  of  zeros  of  h(z)  which  lie  within  the  contour  is  (2n)"^  times  the 
change  in  phase  of  h(z)  as  z  traverses  the  contour.  Now  the  chauige  in 
phase  is 


j.0in  Re**^"*^ 

[arg  [arg  +  [arg  h(z)]j^^i„  +  [arg  h(z)]^^.i^  .  (l.lO) 


As  R-#a>  and  r-»0  ,  the  first  two  terms  vanish  since  h(z)  =  li®  h(z) 

^  z-^0  '  '  z-*<o  '  ' 

=  1  .  Now  a  representation  for  b[  xe"^^^  follows  from  (I.6),  and  so  the 

last  two  terms  of  (l.9)  become 


2n  lim 
r-*0 

R-»as 


arc 


tan 


xe~2x 


l-xe’^^i(2x) 


=  0 


(1.11) 
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Thus  g(z)  has  no  zeros  in  the  complex  plane  cut  along  the  negative  real 
axis. 


It  readily  follows  that  the  path  of  integration  for  (l,3)  may  be  de¬ 
formed  into  the  path  QAP  vd.th  R-*»  and  that  the  value  of  this  integral 
taken  over  the  complete  circuit  APBCDEFQIV  is  zero.  The  integral  over  the 
arcs  PB  and  FQ  tend  to  zero  as  R— ,  and  so  we  have  the  loop  integral 
repre  sentation 


(of) 

^  f  eP^G(p)dp  ,  (1.12) 


(crt-) 

vThere  J  means  that  the  path  of  integration  is  FEDCB  with  r^O  and 
-00 

R-»-  .  We  readily  find  the  Integral  representation 


g(z) 


i  /  xe-^^^^^)dx 

®  l-xe‘^^i(2x)  j  +n^x^e“*^ 


,  Re(z)  >  -2 


(1.13) 


where  Re  stands  for  the  real  part. 

At  this  point,  it  is  helpful  to  give  a  synopsis  of  the  material  in 
the  remaining  sections.  In  Section  II, we  develop  the  usual  power  series  solu¬ 
tion  of  g(z)  for  small  z  ,  and  an  asyugjtotic  ei^ansion  of  g(z)  for  large 
z  .  It  is  the  method  of  constructing  the  power  series  solution  that  motivates 
o\ir  idea  of  using  rational  approximations  to  approximate  G(p)  which  in  turn 
leads  to  an  approximation  of  g(z)  .  Hie  method  of  approximation  including  a 
complete  error  analysis  is  the  subject  of  Section  III.  Numerical  exairples  are 
given  in  Section  IV. 
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II.  THE  CONVEMTIONAL  SOLUTIONS  FOR  SMALL  AND  LARGE  VARIABIE 


To  get  an  ejcpansion  in  ascending  powers  of  z  ,  we  combine  (1.7)  and 
(1.2)  and  formally  develop  G(p)  as  a  series  in  reciprocal  powers  of  p  , 
Termwise  inversion  then  produces  the  desired  result.  Numerous  examples  illus¬ 
trating  this  technique  have  been  given  by  Goldstein  [t].  See  also  Carslaw  and 
Jaeger  [s].  Vfe  find  that 


6(z) 


1  23 

*  k^o  2^kJ 


k-1 

ak  =  (ktl)i  +  ^  (k-r).'ay  . 

r=o 


(2.1) 


Further  numerical  coefficients  are  given  in  the  following  table. 


i- 

^k 

k 

^k 

1 

3 

6 

6431 

2 

11 

7 

62391 

3 

47 

0 

5  24495 

4 

231 

9 

49  60775 

5 

1303 

10 

522  23775 

The  series  (2.1)  converges  quite  slowly.  Indeed,  it  is  clear  from  (l.l3)  and 
a  theorem  in  complex  veu-iable  theory  that  (2.1)  converges  only  if  (z|  <  2  . 
Thus,  if  ^  ^  ^  ^  terms  of  (2.1),  the  result  is  only  correct  to 

about  five  decimals.  Again,  at  z  =  1  ,  we  only  get  accuracy  to  about  2 
decimals . 


To  get  an  asyiiq)totic  expansion  of  g(z)  ,  we  use  (l.4)  to  develop  an 
ascending  series  for  G(p)  ,  and  then  employ  (l.l2).  For  other  exan^iles  of  this 
approach  see  Carslaw  and  Jaeger  [s],  and  Ritchie  and  Sakakura  [s].  We  have 
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and 


4G(p)  =  1  +  pe^(ln  p+e)  -  p^e^  |2  +  e^(ln  pt-e)^ 


^  3^2p 

+  p  e 


1  +  4e^(ln  p4-e)  +  e^^(ln  p<-e)^ 


+  . . . , 


e  =  V  +  In  2 


(2.3) 


Nov 


(04-) 


by  =  -i_  P  e^Pp“-ldp  =  ,  R(u)  a  0  , 


2T7i 


nX‘ 


bQ  =  1  ,  bjj  =  0  ,  n  a  positive  integer  , 


(2.4) 


(cH-) 


■=“  '  air  / 


v=n 


'In 


=  (-)  (n-l): 


'2n 


C3,  =  (:lin-i): 
X" 


3  |t(n)-ln  x|  +  3t'(n)  - 


,  etc.  ,  (2.5) 


where  ^f(z)  is  the  logarithmic  derivative  of  the  gamma  function  r(z)  and 
<|r'(z)  is  the  derivative  of  ♦(z)  .  The  union  of  (2.3)  and  (l.l2)  gives 


g(z) 


I  -  Infktl") 

1  *2  yg/^ 


4(z+2)‘^  (z+4) 

.2 


2(z+6)^ 

+  . . . ,  0  <  z  — .  . 


(z+4)^ 

+2  -  ii  In 


(2.6) 
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Note  that  (2.6)  is  not  very  effective  unless  z  is  quite  large,  say  z  >  6  . 
Also,  the  formula  is  deficient  since  a  useful  estimate  of  the  error  is  not  at 
hand,  although  this  drawback  is  not  too  serious. 

Hie  representations  (2.1)  and  (2.6)  are  therefore  only  useful  for  z 
very  small  or  z  very  large.  There  is  need  for  analytical  representations  to 
cover  the  intermediate  range.  We  shall  show  how  rational  approximations  to 
-Ei(-z)  can  be  effectively  and  economically  used  to  approximate  g(z)  for  z 
ranging  from  the  origin  into  the  region  where  use  of  the  asymptotic  expansion 
is  meaningful. 


III.  SOLUTION  BY  RATIONAL  APPROXIMATIONS 


It  is  known  that  the  exponential  integral  may  be  approximated  by  a 
sequence  of  rational  functions  commonly  known  as  the  Padfe  approximates  or  the 
Gaussian  convergents.  These  approximations  are  related  to  the  asymptotic  ex- 
peinsion  for  -Ei(-z)  ,  see  (l.7),  and  for  the  same  number  of  terms  give  much 
more  accuracy  than  the  asymptotic  expansion.  Indeed,  the  rational  approximates 
can  be  used  effectively  in  regions  \diere  the  asymptotic  expansion  is  useless. 
Thus,  if  the  Pad6  approximates  are  used  to  represent  G(p)  ,  then  an  approxima¬ 
tion  to  g(z)  is  readily  composed,  and  we  should  expect  results  valid  in  a 
much  larger  domain  of  z  values  than  that  derived  using  the  asymptotic  expan¬ 
sions.  This  is  truly  the  case,  and  we  now  take  up  the  procedure  using  rational 
approximates. 

In  [l],  we  showed  that  the  incomplete  gamma  function  which  includes 
the  exponential  integral  as  a  special  case  can  be  approximated  by  a  sequence 
of  rational  approximates  together  with  a  remainder  term.  An  asymptotic  esti¬ 
mate  of  the  error  was  developed  auid  used  to  prove  that  if  the  variable  z  is 
fixed,  jarg  z  |  <  n  ,  then  the  sequence  of  rational  approximates  converge. 

That  is,  the  error  approaches  zero  as  the  order  of  the  approximation  n  is 
increased.  The  asymptotic  estimate  of  the  error  is  qvdte  realistic.  However, 
the  error  formulation  suffered  frcxn  the  fact  that  it  was  not  uniform  in  z  . 
This  deficiency  was  corrected  in  a  recent  study  (]lo3 .  We  sximmarize  below  our 
findings  for  the  exponential  integral.  As  previously  remarked,  similar  results 
are  available  for  the  incomplete  gamma  function  in  the  reference  cited.  We 
have 


-ze^Ei(-z)  =  Vjj(z)  +  Tjj(z)  , 


Vn(^) 


Fn(z) 


Tn(2) 


Un(z) 


f 


(3.1) 


8 


where 


n 


En(z)  =  ^ 
k=o 


=  E  (-)V  E 


(-)’'(y  (r-k)J 
(r+1): 


(3.2) 


and  Fn(^)  2^2  (3.2)  with  k=o  vAence 

Fn(2)  =  lFi(-n;2;-z)  (3.3) 

which  when  suitably  normalized  is  the  Laguerre  polynomial  I^^)(-z)  .  Also 


oo 

Ujj(z)  =  -ze^  (t-z)”f'“"^e"*dt  =  -n.'<((n+l,0;z)  (3.4) 

z 


where  t(a,c;z)  is  a  confluent  hypergeometric  function,  and 


Un(z)  -n.'z’'^'^2f’o(''‘*'l»°‘^2;-z"^)  ,  |z|-*®  ,  |arg  z  |  <  3n/2 


Both  Ejj(z)  and  Fjj(z)  satisfy  the  recurrence  equation 


(3.5) 


(n+2)Fjj+i(z)  =  (z+2n+2)Fn(z)  -  Fn_j_(z) 


(3.6) 


Some  further  useful  relations  are 


zD^+(z+2)D-nJ  Fjj(z)  =0  ,  D  =  d/dz  ,  (3.7) 
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zDFn(z)  =  n  |^Fn(z)-Fn.i(z)J 


(3.8) 


zDEn(z)  =  (z+n+l)Ejj(z)  -  nE^.^Cz)  -  zFn(z)  .  (3.9) 


Let 


z  =  4(n+l)Einh^  a 


(3.10) 


where  a  is  real  and  positive  if  z  is  real  and  positive.  To  be  explicit,  it 
is  convenient  to  write 


sinh^  a 
0. 


peiP  , 

P 

>  0  , 

1  e  [  <  n  , 

3+16  , 

3 

and  6 

real  ,  P  >  0  , 

1  6  1  <  n/2  , 

(3.11) 


so  that 


cosh  f)  = 


^-4p  sin^ 
~2 


} 


sin  5  cosh  P  = 


sin  9/2 


1 

cos  6  sinh  P  =  cos  0/2  .  (3.12) 


Then 


-2rze 


z-2(n+l)(2c+sinh  2a.) 


Tn(z) 


.  Pl(a)  Pp(a)  .3 

J.+  — -  +  — - +0(n  ) 

(r>l)g 


'1-  !!<:> .  ^  .o(„-h 

L  (n+l)2 


Fl(c')  =  (S6  coth^  c)  ^  j^9  coth^a  -6  cot'a^o,  +5^  ,  (3.13) 

Equation  (3.13)  concluded  on  next  page. 
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PgCo)  =  (18432  coth®  c)  j^-135  coth®cr  +180  coth^er.  +558  coth^p  -924  coth^a  +385^, 

(3.13) 


uniformly  in  z  for  n  large  and  |arg  zj  <  n  .  Under  these  conditions  it 
is  clear  that 


lim  Tn(z)  =  0  .  (3.14) 

n-»» 


Thus  the  rational  approximation  cMiverges  every/diere  in  the  con5)lex  plane 
except  when  z  is  on  the  negative  real  axis.  This  is  to  be  expected  since 
the  zeros  of  the  Laguerre  polynomial  lie  on  the  negative  real  axis, 

and  obviously  the  rational  approximates  cannot  be  used  within  sufficiently 
small  neighborhoods  of  these  zeros.  As  a  remark  aside,  in  [lO]  it  is  shown 
that,  if  n  is  fixed,  z  =  -x  ,  and  x  is  sufficiently  large,  then  Vjj(-x) 
may  be  used  to  approximate  xe”*Ei(x)  .  Thus  the  asymptotic  expansion  for 
xe''^Ei(x)  ,  see  (l.7),  and  Vjj(-x)  exhibit  the  same  behavior.  Indeed,  with 
n  fixed,  we  can  show  that 


lim  x^n(.x)  =  0  . 

x-+» 


This  follows  from  an  asymptotic  representation  for  Tjj(-x)  which  is  similar 
in  form  to  (3.13).  For  further  details,  see  [lO^. 

Another  rational  approximation  for  the  exponential  integral  is  given 
by 


-ze^Ei(-z)  =  S^(z)+Bjj(z)  , 


Sn(z) 


_  Mn(^) 
Nn(^) 


Rn(2) 


Wn(z) 

N„(z) 


(3.15) 
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\rtiere 


Mn(z) 


n-1  /  V 

Y. 

k=o  (k+i):(k+i) 


2^2 


/-n+k+1,  1 
2+k,2+k 


n  n 

>  I  (-)V  Z 

k=l  r=k 


(-)^(;)(r-k)i 

Tl 


> 


and  Njj(z)  is  the  2F2  in  (3.16)  vith  k  =  -1  ,  so  that 


Nn(z)  =  iFi(-n;l;-2)  . 


Also 

a> 

Wn(z)  =  ze^  J'  (t-z)’^"^'^e"'^dt  =  nJzt|i(n+l,l;z)  . 
z 

Both  Mji(z)  and  Nn(z)  satisfy  the  recurrence  equation 

(n+l)Nn+3^(z)  =  (z+2n+l)Nj^(z)-nMjj.i(z)  , 


(3.16) 


(3.17) 


(3.13) 


(3.19) 


and  (3.9)  is  valid  if  En(z)  and  Fn(z)  are  replaced  by  Mn(z)  and  Nn(z)  , 
respectively.  Also,  Njj(z)  satisfies  (3.8),  and 

j^zD2+(z+l)D-n  jNn(z)  =  0  .  (3.20) 

A  formula  for  Rn(z)  which  resembles  (3.13)  can  be  given..  However,  it  is 
convenient  to  exhibit  a  foimula  idiich  depicts  the  relation  between 
Tjj(z)  .  Thus 


12 


o 


-Tj^(z)exp 


2cf+ 


tanh  c  +  fcanh  a)(gcosh^c»  +  l)  +o(n”^) 
48(n+l)^  cosh^  a 


X 


+ 


coth  a  ^ 
4(n+l) 

1 

64  cosh^  a 


p 

coth  a  L  1 

Pl(c'.)  - 

coth  o 

32(n+l)^  (n+1)^ 

8 

• 

^oth  o-  2tanh  a+  5tanh^  cr!j 

+0(n-2) 

(3.21) 


V7hence  Rn(z)  enjoys  the  same  convergence  properties  as  Tji(z)  . 

With  these  preliminaries  out  of  the  way.  we  now  pass  on  to  applica¬ 
tion  of  the  rational  approximates  of  Ei(-z)  to  approximate  the  solution  of 
(l.l).  Combining  (l.2)  and  (3.1),  we  can  inrite 


G(p)  =  Gn(p)  - 


h  Tn(2p) 


r2+2pe^Ei(-2p) 


L 


2-Vn(2p) 


Gn(?)  =  i 


_  1  [Fn(gp)-En(gp)] 

[2Fj,(2p)-En(2p)] 


(3.22) 


Thus  G  (p)  is  a  rational  approximation  to  G(p)  .  Let  gn(2)  be  the  inverse 
transform  of  Gn(p)  .  Then  gn(z)  ^  approximation  to  G(z)  ,  and 


G(z)  =  gjj(z)+en(z)  / 

C+lo 

r 

C-ioo 


En(2)  =  -(4^ni)  y’ 


e^PTn(2p)dp 


2+2pe^Ei  (  -2p)  1  2+2pe^Ei  ( -2p)+Tn(  2p)  ] 


(3.23) 


The  numerator  and  denominator  polynomials  of  Gn(p)  follow  from  (3.2).  ,  Clearly 
they  are  easy  to  generate  since  each  satisfies  (3.6)  with  z  replaced  by  2p  . 
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is  readily  decomposed  into  partial  fractions,  and  by  inversion 
is  a  sum  of  exponential  functions.  Thus 


n 


(^)  =  E 

k=l 


-Xjj.z 


“n 


(0)  =  1/4 


(3.24) 


Another  approximation  follows  from  (3.15),  and  for  this,  the  equations  (3.22)  - 
(3.24)  hold  with  obvious  change  of  notation. 

We  now  prove  that  the  approximation  process  is  convergent.  That  is, 
if  z  is  fixed,  jarg  z  j  <  n  ,  then  lim  c,-j(z)  =  0  .  We  have  already  shown 

n-»» 

that  h(p)  ,  see  (l.9),  has  no  zeros  in  the  complex  plane  cut  along  the  nega¬ 
tive  real  axis,  and  obviously  there  are  no  zeros  in  the  right  half  plane.  Let 
c  be  fixed.  If  c  is  sufficiently  large,  the  asymptotic  expansion  of  h(p) 
which  follows  from  (1.6)  -  (l.7)  and  the  asymptotic  form  of  T^(z)  ,  see 
(3.13),  show  that  2h(p)  +  Tjj(2p)  has  no  zeros  on  or  outside  a  semicircle  of 
radius  c  with  center  at  the  origin  which  lies  in  the  right  half  plane. 

Indeed,  since 


2h(p)  1  +  ^  -  _L.  +  0(p"^)  ,  |p|-^”  ,  |arg  p|  <"  3n/2  , 

^  2p^ 


to  prove  convergence,  it  is  sufficient  to  consider  the  integral 


c+ik 

6n(z)  =  -  i  lira  P  Aj^(p)dp  , 
R->a>  y.ik 


An(p)  =  peP(^-^2)-2(n+l)(2o.+sinh  2a)  ^  ^  ^  .  (3.25) 


14 


Let  be  the  path  of  a  quarter  circle  with  center  at  the  origin  and  radius 
c-j^  s  c  which  lies  in  the  upper  right  half  plane.  Let  C2  be  the  imaginary 
axis  extending  from  the  point  icj^  to  iR  .  Then 


=  2Re 


-i  f  ,  I2  =  2Im ,  lim  J  Ajj(p)dp  I 

Cl  [R-^“  Cg  J 


(3.26) 


where  Re  stands  for  the  real  part  and  Im  stands  for  the  imaginary  part. 
We  consider  Ii  first.  Along  Ci  ,  p  =  cie^®  so  that 


n/2 

II  =  2cf  Re  f  exp  I  2ia4-ci(z+2)e^®-2(n+l)(2a+Binh  2a) 
0 

2  *^1®^^  /o 

sinh*^  a  =  -= - r  >  0  Q  ^  n/2  . 

2(n+l) 


de  , 


(3.27) 


We  now  let  n  be  fixed  but  large  with  respect  to  ci  .  !Rien 


cosh  a  = 


cie 


i0 


1  + 


2(n+l) 


Cie 


=  1  + 


ie  ^2^219 


4(n+l)  32(n+l)^ 


+  . . . 


sinh  2a  =  2 


2a  =  2 


®1® 


ie 


2(n+l) 

ciei8 


1  +  _ +  0(n-2) 

4(n+l) 


r 


»ie 


1  -  .  0(n-2) 

12(n+l) 


(3.28) 
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and 


1. 

2(n+l)(2c,+sinli  2a)  =  4  ^2c3^(n+l)J^ 


.ie 


1  +  _!i! _ +  0(n-2) 

12(n+l) 


Since  0  s  ?  s  tt/2  ,  we  have 


(3.29) 


Re 


1. 

2(n+l)(2a+sinh  2o.)  |  £  4|^CQ^(n+l)]  [l  +  0(n~^)J  ,  (3. 


30) 


€uid  as  z  and  are  fixed,  it  is  clear  that  lin  =  0 

n-»® 

For  I2  t  we  let  p  =  iq  along  C2  •  Then 


Ig  =  -2Im 


lim  J  qeiq(2+2)-2(n+l)(2o+sinh 


R-*oo  c 


sinh^  a  =  Si 


iTT/2 


2(n+l) 


(3.51) 


We  suppose  that  0  ^  arg  z  <  rr  so  that  in  (3.3l),  is  always 

bounded  for  C]_  ^  q  ^  ™  .  It  is  convenient  to  split  the  range  of  integration 
into  two  parts.  The  first  from  c^  to  =  4(n+l)  and  the  second  from  R-j^ 

to  R  .  For  convergence,  it  is  therefore  sufficient  to  consider  the  integrals 


I3  =  Re 


Xa.  =  Re 


^g-2(n+l)(2a+sinh  2c,)^^ 


1.^1 


J  ^^dq 


R1 


(3.32) 
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Now  apply  (3.11)  -  (3.13)  irtth  0  =  tt/2  and  p  =  c[/2(n+l)  .  Then 


2(n+l)Re(sinh  2b)  =  2[(n+l)qj‘ 


(l+p2)  -p 


e-2(n+l)Re(2a) 


(l+p2)2+p+2i"p2 


(l+p2)2+p 


i.n-2n-2 
2 


For  the  evaluation  of  I3  ,  0  < 


2(n+l) 


s  p  ^  2 


In  thifj  case 


.-2(n+l)Re(2a)  < 


1  + 

^1 

2 

.  2(n+l)_ 

-8(n-t-l) 


1  + 


2(n+l) 


j2(n4-l)J^  (.4)[2ci(n+l)] 


Now  11m  (l+x"^)  =  ,  and  so 

X-*a> 


lim  I3  <  ^  lim  ri6(n+l)^-c^  1  e"^  [2ci(n+l)]^  ^  Q 

n-»<»  ^ 


(5.33) 


(3.34) 


(3.35) 


(3,36) 
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In  ,  let  and  put  q  =  2(n+l)t  .  Then 


=  4(n+l)%e 


^g-2(n+l)(2cf-+sinh  acr,)^^ 


sinh^  a  =  . 

(3.37) 


Thus  (3.33)  and  (3.34)  apply  with  p  =  t  and  2  s  t  ^  ®  .  In  this  case 

g-2(n+l)Re(2p.)  ^  ^^^j-2n-2  _  (3.38) 

Also  2(n+l)Re(sinh  2o.)  ^0  for  all  t  ,  2  s  t  s  «■  .  Hence 

I .  <  (n-H)^  r  dt  ^  (n+1)^ 

22n*2  ^  ^2n+l  n.2^'^'^2 

Um  I4  »  0  ,  (3,39) 

n-*«* 


and  the  proof  of  convergence  is  complete.  A  similar  argument  shows  that  the 
approximation  technique  based  on  (3.15)  also  converges.  In  this  connection, 
see  (3.21). 

We  have  attempted  to  develop  a  useful  and  realistic  approximation  for 
ejj(z)  so  that  one  could  determine,  a  priori,  for  a  given  z  the  value  of  n 
required  to  achieve  a  given  level  of  accToracy.  Some  preliminary  results  in 
this  direction  are  available,  but  the  analysis  is  far  from  complete,  and  we 
defer  further  comments  for  a  later  report.  It  is  pertinent  to  remark,  however, 
that  numerical  values  determined  by  the  approximation  technique  can  be  used  to 
give  pragmatic  estimates  of  the  error.  This  and  much  more  are  brought  out  in 
the  follo^d.ng  section. 


t 
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In  Table  1,  ve  tabulate  J  f^e'^dt  ,  its  rational  approximation 

z 

(see  (3.1)),  the  exact  error,  and  the  approximate  error  according  to  (3.13). 

The  data  are  developed  for  z  =  and  n  =  4  .  Note  that  the  error  formula 

for  Tjj(z)  is  quite  remarkable  for  it  is  very  realistic  even  for  small  n  . 


TABLE  1 


z'le"Zv^(z) 


0 

"1/4 

rr/2 

3t:/4 


4.89005*10-2 

.3 . 95846 *  10 "2-8 . 22921i * 10"2 
-4. 22981*10"^3. 461671*  10-2 
-2.16947  +0.317771 


4.89190*10-2 

-3 . 95652 • 10-2-8 . 224141  *  10-2 
-4.23980*10-1+3 . 436821 • 10-2 
-2.12827  +0.364621 


z'^e"'*T4(z) 

z-le-’^T4(z) 

Exact 

Approximate 

0 

-0.185*10-^ 

-0.184*10-'J^ 

n/4 

-0.194*10-^-0.5071*10-^ 

-0.194*10-^-0.507i*10-4 

tt/2 

0.999* lO-^+O. 246i * 10"^ 

0.996*10-3+0. 249i *  10-3 

3n/4 

-0.412*10-1-0. 468i *10-1 

-0 . 412* 10-1-0. 473i • lO"! 

In  Table  II,  we  record  the  values  of  aj^  and  for  n  =  3(l)6  , 

and  in  Table  III  we  compare  gjj(z)  with  the  "exact"  g(z)  for  several  values 
of  z.  These  data  are  based  on  the  representation  (3.1).  Tables  17  and  V  are 

the  same  as  I  and  II,  respectively,  save  that  the  approximations  are  based  on 

(3.15).  Here  we  use  the  notation  ajj*,  and  gjj*  ,  which  parallels  that 

of  (3.24).  The  "exact"  values  for  z  s  5  are  taken  from  a  study  by  Fox  and 
Goodwin  [ll],  \dio  use  (l.l)  in  a  study  of  numerical  solutions  to  integral  equa¬ 
tions.  The  value  of  g(z)  for  z  =  6  was  found  from  (1.13)  by  ntunerical 

integration.  It  appears  that  the  value  of  g(z)  for  z  =  5  is  a  unit  too 
large  in  the  last  place.  In  Tables  III  and  V,  we  also  give  values  of  8(2) 
for  z  =  4(1)6  -vrtiich  are  derived  from  the  n  =  4,5  and  6  approximates  using 
Aitken's  S^-process. 
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TABI£  II 


II 

n  =  4 

k 

% 

^k 

k 

ak 

^k 

1 

0.09609  77 

-0.71098 

97 

1 

0.05937  40 

-0.53891 

09 

2 

0.14655  50 

-1.89746 

62 

2 

0.15682  37 

-1.56937 

11 

3 

0  -  0073  4  73 

-3.89154  41 

3 

0.03138  23 

-2.91402 

74 

4 

0.00042  00 

-5.47769 

05 

n  =  5 

n  =  6 

k 

^k 

^k 

k 

ak 

^k 

1 

0.03979  98 

-0.42941 

92 

1 

0.02843  62 

-0.35499 

71 

2 

0.14295  32 

-1.33416 

55 

■'2 

0.11665  90 

-1.15071 

93 

3 

0  -  06412  69 

-2.40159 

29 

3 

0.09292  31 

-2.07864 

32 

4 

0.00309  96 

-4.20473 

33 

4 

0.00973  90 

-3.47492 

24 

5 

0.00002  05 

-7.13008 

91 

5 

0.00024  17 

-5.61773 

46 

6 

0.00000  10 

-8.82298 

34 

z 

83(2) 

64(z) 

TABUE  III 

65(2) 

g6(=^) 

6(z) 

6(z) 

0 

0.25000 

0.25000 

0.25000 

0.25000 

0.25000 

0.5 

0.  12515 

0.12515 

0.12515 

0.12515 

0.12515 

1.0 

0.06932 

0.06941 

0.06941 

0.06941 

0.06942 

1.5 

0-04161 

0.04193 

0.04197 

0.04197 

0.04197 

2.0 

0  .  02648 

0.02716 

0.02730 

0.02732 

0.02733 

3.0 

0.01138 

0.01322 

0.01363 

0.01374 

0.01378 

0.01378 

4.0 

0.00567 

0.00718 

0.00783 

0.00809 

0.00826 

0.00821 

5.0 

0  .  00267 

0.00407 

0.00483 

0.00520 

0.00555 

0.00547 

6.0 

0.  00135 

0.00235 

0.00307 

0.00350 

0.00414 

0.00391 
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TABIE  IV 


n  =  3 

II 

* 

\  * 

* 

\  * 

k 

^k 

k 

^k 

^k 

1 

0.03668  99 

-0.31385  -93 

1 

0.02025  70 

-0.22504  61 

2 

0.18750  00 

-1.50000  00 

2 

0.15173  03 

-1.18873  16 

3 

0.02581  01 

-3.18614  07 

3  = 

0.07630  34 

-2.38574  42 

4 

0.00170  92 

-4.70047  81 

n  =  5 

n  =  6 

1  * 

■V  * 

k 

^k 

k 

^k 

1 

0.01279  44 

-0.17429  89 

1 

0.00880  86 

-0.14178  99 

2 

0.11138  06 

-0.96597  23 

2 

0.08145  38 

-0.80353  59 

3 

0.11584  45 

-1.98033  70 

3 

0.13290  86 

-1.71581  14 

4 

0.00989  00 

-3.55882  82 

4 

0.02597  15 

-2.92720  02 

5 

0.00009  05 

-6.32056  36 

5 

0.00085  32 

-4.92022  14 

6 

0.00000  43 

-7.99144  12 

z 

e3*(-) 

64*(z) 

TABUS  V 

%*(-) 

g6*(z) 

6*(z) 

6(z) 

0 

0.25000 

0.25000 

0.25000 

0.25000 

0.25000 

0.5 

0.12518 

0.12515 

0.12515 

0.12515 

0.12515 

1.0 

0.06971 

0.06943 

0.06941 

0.06941 

0.06942 

1.5 

0.04289 

0.04210 

0.04197 

0.04197 

0.04197 

2.0 

0.02896 

0.02764 

0.02738 

0.02733 

0.02733 

3.0 

0.01639 

0.01466 

0.01403 

0.01384 

0.01376 

0.0137G 

4.0 

0.01092 

0.00955 

0.00875 

0.00641 

0.00813 

0.00821 

5.0 

0.00774 

0.00697 

0.00625 

0.00582 

0.00524 

0.00547 

6.0 

0.00560 

0.00537 

0.00484 

0.00442 

0.00282 

0.00391 
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As  previously  remarked,  our  formulation  of  the  error  in  the  approxi¬ 
mation  process  is  not  in  a  convenient  fom  for  practical  use.  Since  the  ap¬ 
proximate  solutions  converge  and  are  easy  to  generate,  it  is  suggested  that 
the  numbers  themselves  be  employed  to  apprise  the  accuracy  of  the  process. 

Thus  if  two  or  more  successive  convergents  are  evaluated,  one  can  accept  the 
common  digits  as  correct.  Further,  as  shown  in  the  above  tables,  Aitken's 
6^-process  is  valuable  provided  n  ,  the  order  of  the  approximation,  is  suf¬ 
ficiently  large.  Note  that  in  Tables  II  and  IV,  the  62_process  gives  improved 
values  except  for  z  =  6  in  Table  IV. 

Another  technique  useful  to  assess  the  error  is  to  compare  the  ap¬ 
proximate  entries  for  large  z  with  values  deduced  from  the  asymptotic  expan¬ 
sion.  If  we  use  all  the  terms  given  in  (2.6),  then  for  z  =  5  and  6  ,  we  get 
0.00509  and  0.00374,  respectively.  In  our  present  example,  the  asymptotic 
representation  is  deficient  since  an  error  estimate  is  not  available,  and  it 
is  clear  from  the  calculations  that  the  asymptotic  expansion  is  not  suffi¬ 
ciently  reliable  unless  z  is  considerably  larger  than  6.  Of  course,  had  we 
computed  approximates  for  larger  n  ,  then  we  could  take  z  larger  whence 
values  deduced  from  the  asymptotic  representation  would  be  more  informative. 

In  our  study,  a  further  procedure  to  get  improved  values  is  to 
average  the  corresponding  entries  in  Tables  III  and  V.  This  works  because  in 
Table  III  the  values  are  converging  from  below  while  in  Table  V,  the  values 
are  converging  from  above.  Thus  for  z  =  5  and  6  ,  we  get  the  improved  values 
0.00551  and  0.00396,  respectively.  Of  course,  in  practice  it  may  not  always 
be  possible  to  have  rational  approximations  so  that  one  can  get  convergence 
from  both  sides.  In  any  event,  it  seems  that  using  only  one  sequence  of 
rational  approximates  \7ith  more  entries  for  n  would  be  more  informative. 


22 


REFERENCES 


1.  Luke,  Y.  L. ,  "The  Padfe  Table  and  the  T-Method,"  Jour.  Math.  Phys.,  37, 

110-127  (1958). 

2.  Luke,  Y.  L. ,  "On  Economic  Representations  of  Transcendental  Functions," 

Jour.  Math.  Phys.,  38,  279-294  (i960). 

3.  Luke,  Y.  L. ,  "On  the  Approximate  Inversion  of  Sore  Laplace  Transforms," 

Aeronautical  Research  Laboratory  Report  ARL  20  (May,  1961).  This  is  also 
to  be  published  in  the  Proceedings  of  the  Fourth  U.  S.  National  Congress 
of  Applied  Mechanics. 

4.  Luke,  Y.  L. ,  "Approximate  Inversion  of  a  Class  of  Laplace  Transforms  with 

Application  to  Supersonic  Flow  Problems,"  Midwest  Research  Institute 
Report,  July  1961. 

5.  Friedlander,  F.  G. ,  "The  Reflexion  of  Sound  Pulses  by  Convex  Parabolic 

Reflectors,"  Proc.  Cambridge  Philos.  Soc.,  134-149  (1941). 

6.  Erdfelyi,  A.,  Magnus,  V.,  Oberhettinger,  F. ,  and  Tricani,  F.  G. ,  Higher 

Transcendental  Functions,  v.  1,  McGraw-Hill  (1953). 

7.  Goldstein,  S.,  "Two-Dimensional  Diffusion  Problems  with  Circular  Syumetry, " 

Proc.  London  Math.  Soc.  (2)  34,  51-88  (1932). 

8.  Carslaw,  H.  S.,  and  Jeieger,  J.  C.,  "Conduction  of  Heat  in  Solids,"  Oxford 

pp.  330-342  (1959). 

9.  Ritchie,  R.  H. ,  and  Sakakura,  A.  Y. ,  "Asymptotic  Expansions  of  Solutions  of 

the  Heat  Conduction  Equation  in  Internally  Bounded  Cylindrical  Geometry," 
Jour.  Appl.  Phys.,  27,  1453-1457  (1956). 

10.  Luke,  Y.  L. ,  "Rational  Approximations  for  the  Incomplete  Gamna  Function," 

Midwest  Research  Institute  Report,  July  1962. 

11.  Fox,  L. ,  and  Goodwin,  E.  T. ,  "The  Numerical  Solution  of  Non-Singular  Linear 

Integral  Equations,"  Philos.  Trans.  Roy.  Soc.  London.  A  245,  501-534 
(1953). 


23 


